Finite-density corrections to the Unitary Fermi gas: A lattice perspective from Dynamical 

Mean-Field Theory 
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We investigate the approach to the universal regime of the dilute unitary Fermi gas as the density is reduced 
to zero in a lattice model. To this end we study the chemical potential, superfluid order parameter and internal 
energy of the attractive Hubbard model in three different lattices with densities of states (DOS) which share the 
same low-energy behavior of fermions in three-dimensional free space: a cubic lattice, a "Bethe lattice" with 
a semicircular DOS, and a "lattice gas" with parabolic dispersion and a sharp energy cut-off that ensures the 
normalization of the DOS. The model is solved using Dynamical Mean-Field Theory, that treats directly the 
thermodynamic limit and arbitrarily low densities, eliminating finite-size effects. At densities of the order of 
one fermion per site the lattice and its specific form dominate the results. The evolution to the low-density limit 
is smooth and it does not allow to define an unambiguous low-density regime. Such finite-density effects are 
significantly reduced using the lattice gas, and they are maximal for the three-dimensional cubic lattice. Even 
though dynamical mean-field theory is bound to reduce to the more standard static mean field in the limit of 
zero density due to the local nature of the self-energy and of the vertex functions, it compares well with accurate 
Monte Carlo simulations down to the lowest densities accessible to the latter. 

PACS numbers: 71.10.Fd, 03.75.Ss, 05.30.Fk, 71.10.-w 



I. INTRODUCTION 

The ever increasing ability to manipulate and control ultra- 
cold Fermi gases allows to experimentally realize in an accu- 
rate and controlled way physical conditions which can only 
be approximately realized in solid state, like the crossover 
from weak-coupling Bardeen-Cooper-Schrieffer(BCS) super- 
fluidity to strong-coupling Bose-Einstein condensation (BEC) 
of preformed bosonic pairs, which occurs in two-component 
fermionic systems as a function of the coupling strength'. 

Actual realizations of degenerate ultracold Fermi gases are 
dilute, i.e., the mean interparticle distance is much larger than 
the range i?o of the interatomic potential. In terms of the 
Fermi momentum kp = (37r^n) 3 the diluteness condition is 
equivalent to kpRo <C 1. For such dilute systems the effect of 
the interaction can be parameterized through the s-wave scat- 
tering length Os only, regardless the details of the potential, 
provided its short-range character. This is particularly intrigu- 
ing because Fano-Feshbach resonances permit to control the 
scattering length as by simply tuning a magnetic field. One 
can thus move from negative values on the BCS side to pos- 
itive values on the BEC side, passing through the resonance 
point, where diverges. This latter situation is usually re- 
ferred to as unitary limit, and it displays an extra universality 
because the divergence of leaves us with a single length 
scale cx n~^/^ and a single energy scale, the Fermi energy^. 
For example, the chemical potential and the superfluid energy 
gap will depend on the density only through the Fermi en- 
ergy Ep, i.e. fi,A(x Ep though with non trivial coefficients. 
The evaluation of these coefficients defies simple analytical 
treatments, due to the strongly coiTelated nature of the prob- 
lem, and to the lack of obvious expansion parameters. In this 
light numerical methods, like Quantum Monte Carlo (QMC) 
simulation s^i'^i^'^'^-^-^' "' of finite systems have an extremely im- 
portant role, both to obtain direct estimates of the observables. 



and to guide the choice of the relevant classes of diagrams in 
perturbative expansions. The use of QMC overcomes the lim- 
itations of perturbative methods, but QMC simulations still 
suffer from finite-size effects. Moreover they require a spe- 
cific choice of the interaction potential, which may be relevant 
if the density is not small enough to assure the realization of 
the dilute limit. 

An alternative to a direct simulation of the dilute Fermi gas 
is to consider the low-density limit of three-dimensional lat- 
tice models with local interaction, like the attractive Hubbard 
model. Obviously the lattice introduces a new length scale 
associated with the inter-site spacing. Therefore the univer- 
sal behavior can only be recovered in the zero density limit, 
where the infinite inter-particle distance makes the presence 
of the lattice iiTelevant. This approach has been used in the 
QMC simulations of Refs. 5.7 8 9, and it has the advantage of 
an intrinsic regularization given by the lattice spacing, which 
from the diagrammatic point of view introduces an ultravio- 
let cutoff. Nonetheless, even if QMC results for the attractive 
Hubbard model can be regarded as essentially exact for finite 
systems, it should be kept in mind that actual simulations are 
limited to a finite number of sites and the density can not be 
arbitrarily reduced. Hence they are plagued by both intrin- 
sic finite-size effects and by finite-density effects, calling for 
careful extrapolations to the thermodynamic and zero-density 
limits. 

In this paper we follow the lattice route to the dilute limit 
using a theoretical approach which is explicitly built in the 
thermodynamic limit and it has no intrinsic limitations to treat 
arbitrarily small densities, the Dynamical Mean-Field Theory 
(DMFT)". DMFT is a quantum, and more accurate, version 
of standard mean-field theories, in which the solution in the 
thermodynamic limit is obtained neglecting spatial coiTela- 
tions. This approach has been widely used in the context of 
solid-state physics'', and it proved accurate for most three- 
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dimensional correlated solids, and more recently it has been 
proposed to access the properties of the dilute Fermi gas in 
Ref. [13, where the system was studied in a three-dimensional 
cubic lattice. Since our focus is to understand the relevance 
of non-universal corrections to the universal regime, here we 
extend the analysis of Ref. [13 considering three different lat- 
tices (which will have different finite-density corrections), and 
much smaller densities. 

Our analysis shows that the extrapolation to zero density 
requires extreme care. Even if the three lattices have the same 
low-energy density of states (which yields the same zero den- 
sity limit), they provide significantly different results (i.e., 
they are not in a universal regime) down to densities of the 
order of 0.01 fermions per site. This is result is an important 
warning because present QMC calculations are limited to den- 
sities of this order of magnitude (to our knowledge the lowest 
density available for the present problem is n = 0.05 in Ref. 
i). 

The weakest dependence on density is found for a model, 
the so-called lattice gas, in which the density of states (DOS) 
coincides with that of fermions in free space up to a cut-off 
after which it vanishes. The results for the cubic lattice and 
for a lattice with a semicircular DOS are influenced to a larger 
extent by the large-density behavior Only for extremely small 
densities the three lattices provide the same result. Unfortu- 
nately we find that in the same small-density limit DMFT re- 
duces to static mean-field (MP), since the method is unable to 
reproduce the divergence of beyond MR However we will 
see that the limitations of DMFT appear only at very small 
densities, while DMFT introduces important corrections to the 
static mean-field down to densities well below n ~ 0.05. This 
observation also explains why in Ref. 12 DMFT was found to 
extrapolate to a different value with respect to static MF. In- 
deed our analysis shows that in this reference the zero-density 
limit has been extrapolated from densities larger than those 
for which DMFT approaches static mean-field. 

The paper is organized as follows: In Sec. II we intro- 
duce the attractive Hubbard model and the approach to the 
zero-density limit. Sec. Ill presents results of static mean- 
field; Sec. IV is devoted to the implementation of Dynamical 
Mean-Field Theory for the attractive Hubbard model; Sec. V 
presents the DMFT results and Sec. VI contains the final re- 
marks. 



II. ATTRACTIVE HUBBARD MODEL APPROACH TO 
THE UNITARY LIMIT 

In this section we discuss how the three-dimensional attrac- 
tive Hubbard model can be used to describe the properties of 
a two-component dilute gas. The Hamiltonian of the model 
reads 

Z!(^fc~^)'^fc.<T + ^Z!'^'.T"^i. (1) 

fc,o- » 

where Uka = c\^Ckcj (riia = c\^Cia) is the number operator 
for Fermions of spin a in momentum k (site i). [/ < is the 



strength of the local pairing interaction between fermions with 
different flavor, and is the free dispersion in a chosen lattice 
with lattice spacing I. 

The attractive Hubbard model has been studied in differ- 
ent frameworks as a paradigm for lattice superconductors and 
superfluidsji^ Here we use this lattice model merely as a sys- 
tematic way to approach the universal regime of an interacting 
gas. Therefore we focus on the small-density regime, where 
the effect of the lattice is bound to disappear because the aver- 
age interparticle distance will ultimately become much larger 
than the lattice spacing. Nonetheless, at every finite density 
the details of the lattice are in principle relevant. 

In our theoretical approach the lattice under consideration 
and its dispersion e^: enter the calculation only through the 
density of states. For a system of non-interacting fermions 
with mass m in three-dimensional continuum space the DOS 
per unit volume is (assuming h — 1 everywhere) Dfreei^) — 

'^^^^ \/e. In order to obtain the correct zero density limit, 
we choose three lattices whose DOS's per lattice site D lattice 
share the same low-density behavior of Dj^ee- This implies 
that Diattice{f-) — Df reei^) for Small energies, where I is 
the lattice spacing that we will use as the unit length. In the 
following we drop the index "lattice". 
Namely we use: 

• A semicircular DOS 

where D is the half bandwidth, that we will use as 
energy unit in the rest of the manuscript. This DOS 
does not correspond to any obvious dispersion but, as 
we shall see, it has some practical advantages. For 
small energies close to the bottom of the band we 
have I?sc(e) ~ \/e, which corresponds to Df^ee 
if we suitably choose the lattice effective mass m = 
{4Tr)^/D. 

• The so-called lattice gas, whose DOS coincides with 
the free DOS up to a cutoff A = 1.40539 .. .D above 
which it vanishes. This value of A ensures the normal- 
ization of the DOS per lattice site. This model can be 
seen as an approximate version of a system in which 
only the momenta included in the Brillouin zone of the 
lattice are allowed, and for each of these momenta the 
dispersion coincides with the parabolic dispersion of 
free fermions electrons. Indeed implementing this pro- 
cedure yields a similar DOS with a smoother energy 
cutoff. This model has been studied in Ref. 151161191, where 
a different value of the cutoff was used. 

• A three-dimensional lattice DOS for nearest-neighbor 
hopping on a cubic lattice with dispersion — 
—2t{coskx + cosky + coskz) with t = l/2m = 
^^^2(47r)^/^. This choice has been used in Ref. [T2I and 

In order to obtain a first insight on the properties of the 
different lattices, in Fig. [T]we compare the low-energy part of 
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FIG. 1: (Color Online) Comparison between the three different DOS's with 
parameters such that the low-energy behavior coincides with that of free elec- 
trons, as discussed in the text: semicircular DOS (dot-dashed red Hne), cubic 
three-dimensional lattice (dashed green line), lattice gas (solid blue line). For 
reference the DOS for free fermions in three-dimensional space is shown as 
a dotted black line. Energies are expressed in units of the half-bandwidth D 
of the semicircular DOS. 



our DOS's (Obviously in this region the lattice gas coincides 
with the target DOS for fermions is three dimensions). It is 
interesting to observe that the semicircular DOS follows very 
closely the free 3-d DOS, and it is therefore likely to allow for 
a good convergence to the dilute limit n 0, while the cubic 
lattice rapidly departs from the reference. 



As far as the interaction strength U is concerned, in order 
to achieve the universal behavior in the dilute limit, U has to 
be chosen as the unitary value, for which the scattering length 
diverges. 

For a Hubbard potential as is related to the coupling con- 
stant U through the relation** 



where 



as{U) 



n(o,o) 



1 



47rft2 [/-I +n(0,0) 



2e 



(3) 



(4) 



For the semicircular DOS one finds Uc ~ —D, while for the 
lattice gas Uc = —1.1624 . . .D and for the cubic lattice Uc = 
-0. 73224,... D. 



III. A FIRST INSIGHT FROM STATIC MEAN-FIELD 

In order to obtain a first feeUng about the role of finite- 
density corrections to the dilute limit of lattice models and to 
understand the role of the precise choice of the lattice, we start 
from a simpler approach with respect to DMFT, i.e. the static 
MF method. This approach shares the main advantages of 
DMFT with respect to exact numerical methods, being im- 
plemented in the thermodynamic limit without any restric- 
tion in terms of attainable densities. Of course the method 
is less accurate than DMFT, which introduces exactly local 



quantum fluctuations, but it can be solved at very low compu- 
tational cost at every density. Finally, the zero-density limit 
of the method is the well-known Leggett's MF theory for the 
BCS-BEC crossover in dilute gases'"*, which we can use as a 
benchmark to quantify finite-density corrections and the de- 
pendence on the actual lattice. 

The static MF approximation consists in decoupling the 
attractive interaction term both in the normal and in the 
anomalous s-wave Cooper channels, determining the value 
of the pairing amplitude and of the chemical potential self- 
consistently. Then we introduce the scattering length using 
Eq. (O. Finally we rescale the relevant energies by the non 
interacting Fermi energy Ep of particles in 3d continuum 
space. The MF equations for the rescaled chemical potential 
fl = fi/Ep and gap parameter A = A/Ep are: 



dx D{x) 



{x - /I') 



(x - ti'Y + A2 



=- 



kpQs 



- dxD{x) 

TT 



{x - fl')^ + A 



(5) 



(6) 



where the left-hand side of Eq. ^ vanishes because we take 
the unitary limit Os ^ oo. D{x) = Tyr^TET)^ is the 
lattice DOS, m' = /i 



is the renormalized chemical po- 
tential including the Hartree contribution (which vanishes in 
the dilute Umit, so it does not appear in Leggett's theory) and 

X = e/ Ep. 

The equations correctly reduce to Leggett's results in the 
zero-density limit, where both the gap and the chemical po- 
tential are proportional Xo Ep. At finite density, corrections to 
the limiting value are introduced both by the different density 
of states with respect to the free DOS used by Leggett, and 
by the Hartree shift of the chemical potential. The numerical 
solution of the MF equations ( I5l6l l for the reduced gap and 
chemical potential is reported in Fig. ^ for the three dif- 
ferent lattices under considerations. In order to highlight the 
effect of the lattice we plot the results from the dilute limit 
71 = to the so-called half-filling {n = 1) situation in which 
we have one fermion per lattice site. The most evident result 
is the slowness of the evolution from n — 1, in which the 
effect of the lattice is maximum and the three models are ob- 
viously different, to n = where the models have to become 
equivalent. The results appear to be strongly influenced by the 
n = 1 value down to very low density, and the smoothness of 
the curves makes it questionable to define a reasonable low- 
density region from the data. As a consequence of the persis- 
tence of the lattice effects, the finite-density behavior remains 
strongly dependent on the choice of the lattice. 

This effect is apparently much stronger on the chemical 
potential, while the superfluid order parameter depends in a 
weaker way on the choice of the lattice. This is most likely 
due to the fact that the unitary regime is physically closer to 
the Bose-Einstein regime, in which the interaction strength 
controls the order parameter and the details of the bandstruc- 
ture are less important, than to the BCS regime, where the 
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in a simulation on a finite lattice, where only a finite number 
of densities can be considered. In other words, the simple ob- 
servation of a zero derivative of the observables as a function 
of density does not guarantee that the system is in the dilute, 
and authentically universal, regime, but rather it can suggest 
an artificial, or fake universality. 



IV. DYNAMICAL MEAN-FIELD THEORY OF THE 
ATTRACTIVE HUBBARD MODEL 

A. The General Formalism and the Model 



FIG. 2: (Color Online) Static mean-field results for the reseated chemical 
potential Ep and superconducting order parameter A/Ep as a function 
of density n. The dashed blue line is for the cubic lattice, the red solid line 
for the semicircular DOS and the dot-dashed green line is for the lattice gas. 
Horizontal arrows mai'k Leggett's results in the zero-density limit. 



DOS at the Fermi level has a major effect on the order pa- 
rameter. Nonetheless, the dependence on density is still very 
strong for all the considered lattices. 

In the case of the chemical potential (and of the renormal- 
ized internal energy, that we will show later in comparison 
with DMFT results), it is evident that the cubic lattice in- 
troduces significantly larger finite-density corrections at any 
density, while the two other models can give rise to smaller 
deviations. As mentioned above, this is essentially due to the 
peculiar value of /2 at n = 1 for the cubic lattice, ultimately 
arising from the peculiar behavior of the DOS. Finally, all the 
curves approach the zero-density limit from below and with 
a very large slope. This is a first indication that any attempt 
to extract informations about the dilute Fermi gas from a lat- 
tice perspective needs a very careful extrapolation to the zero- 
density limit (besides the thermodynamic limit), and that the 
choice of the non-interacting density of states can be crucial 
to minimize these effects. 

Another observation is in order The behavior as a function 
of density of the rescaled quantities derives from a rather in- 
volved interplay between the "large-density" effects (in which 
the precise choice of the lattice is very important) and the ten- 
dency towards the universal unitary limit. As a consequence, 
it is not possible to define a well-defined density scale under 
which the finite-density system already displays the univer- 
sal physics. We can have regimes of density in which some 
quantities appear rather flat, or they have minima. In these re- 
gions or points one has that the derivative of the renormalized 
quantities as a function of the density is zero. For the sake 
of definiteness, we have that, for the semicircular DOS one 
has dfl/dn — 0, at n ~ 0.1. If no data were available for 
smaller densities one could be tempted to interpret the van- 
ishing derivative as the beginning of a regime in which ji is 
independent on n, i.e., a universal regime in which the de- 
pendence on density is lost. Obviously, our MF data show 
that this is not the case, and the vanishing derivative is only 
associated to a minimum after which the dilute limit is ap- 
proached. We observe that this could not be equally evident 



We briefly introduce DMFT and its application to attrac- 
tive models and superfluidity. Previous DMFT studies of this 
model have been so far mostly devoted to the high-density 
regime. In Refs fT5lfT6l the normal phase has been studied 
by excluding superfluidity, and a pairing transition between 
a normal metal and a pseudogapped state of preformed pairs 
has been found. Various properties of the superfluid state have 
been studied with the same approach, highlighting the ability 
of DMFT to properly describe the BCS-BEC crossover with- 
out bias towards one of the two limit o'^'^^i'^i^" . The same ap- 
proach has also been used to study two-component fermionic 
mixtures with mass unbalance^' and density unbalance. 

The low-density regime of the attractive Hubbard model 
has been recently studied within DMFT in order to access 
the unitary regime in the normal state^^ and in the superfluid 
phase'^. In both papers, the cubic lattice has been used, and 
some limitations at low-densities have been reported. Here we 
extend the analysis of Ref. considering different lattices, 
and we improve the numerical accuracy at low density, which 
will allow us to draw important conclusions about the ability 
of DMFT to describe the low-density limit. 

DMFT maps a quantum lattice model onto a local prob- 
lem, which can be represented through an "impurity model"' -, 
i.e., a model in which a single interacting site is embed- 
ded in a non-interacting medium. In our case ([T]l is mapped 
onto an impurity model with attractive coupling, and the non- 
interacting bath is superfluid. Namely 

Ham = ^ [e; 4^ + Vi {c\^d„ + h.c.) 

■Hioc (7) 



La 



where Tiioc — —Unoinoi — fino is the on-site term, and 
the chemical potential /i controls the density. For all the 
different lattices, /i = corresponds to zero density in the 
non-interacting case. From the impurity model we com- 
pute the normal and anomalous Green's functions, G(r) — 
-(Tc|(r)c|(0)) and F{t) = -(Tc|(r)q(0)), and the cor- 
responding normal and anomalous (superfluid) components of 
the self-energy E and S. 

The correspondence between the effective local model 
(|7]l and the original lattice model is guaranteed by a self- 
consistency condition analogous to the Curie- Weiss equation 
for the Ising model. The self-consistency condition can be 
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obtained by requiring that the impurity Green's function coin- 
cides with the local component of the lattice Green's function, 
namely 



G{iuJn) = Glatt{iUJn) = dc D{e) 



z — e 



where z = iujn + M ~ '^{i'-^n)- 

We underline that the original lattice enters only through 
the non-interacting DOS, and that for the semicircular DOS 
dU this equation is particularly simplified, since it does not 
require the numerical integration on the energy e". 

It is important to notice that DMFT can be also seen as an 
approximation of the full Luttinger-Ward functional in which 
only the local Green's function is considered". This means 
that DMFT is a variational method. Since the static MP, in 
which the self-energy is local and it has no dependence on fre- 
quency, is a subcase of DMFT in which the self-energy loses 
the dynamical character, the variational principle implies that 
DMFT is always an improvement over static MF, unless the 
two methods become identicaUi^i. 



B. Exact Diagonalization Solution 

Despite the simplifications introduced in DMFT, the effec- 
tive local model O can not be solved exactly by analytical 
methods and numerical solutions are necessary. Yet, the re- 
quired computational effort is enormously lighter than for the 
original lattice model, and, more importantly, the method is 
defined in the thermodynamic limit, so that no finite-size cor- 
rections need to be considered, as opposed to full numerical 
solution of lattice models. Analogously, we are in the grand- 
canonical ensemble, where we can tune the chemical potential 
and the density can assume arbitrarily small values. 

Here we adopt Exact Diagonalization (ED)-^, which re- 
quires to truncate the sum in Eq. (|7J to a finite number of lev- 
els Ng. In practical implementations Ng is necessarily small, 
but it has been shown that Ng smaller than 10 provides accu- 
rate results for thermodynamic properties. Such a discretiza- 
tion introduces a new step in the iterative procedure. Namely, 
after a new bath is obtained by means of the self-consistency 
equation, it has to be represented into a discrete form. This 
representation is usually obtained by fitting the result of the 
self-consistency to a discrete system. The details of the fit 
have to be chosen with care, especially in "delicate" regions 
such as the low-density regime we are interested in. 

As customary, we perform the fit on the imaginary- 
frequency axis (even if we work at T = 0), where the Green's 
functions are smooth. This requires the definition of a fic- 
tive inverse temperature f3 which defines a Matsubara grid. A 
large (3 is required to investigate the low-frequency behavior 
because of the small energy scales. Here we used values of 
(3 up to 6000 for the smallest densities. All the results we 
present are converged as a function of (3. Other aspects of the 
fit are discussed, e.g., in Ref. I26ll27l 



ED has been used in the low-density limit of the attrac- 
tive Hubbard model in Ref. [l^ where full diagonalization 
of the Hamiltonian matrix has been used, limiting the study to 
A^s = 6 (5 levels in the bath). Here we use the Lanczos al- 
gorithm, which allows to compute the groundstate for a larger 
Ng. We analyzed systematically the behavior of the results for 
Ns = 6, 7, 8, 9, and we found that the results for static observ- 
ables are essentially converged as a function of the number of 
bath levels for Ng = 8, which is the value that we will use 
throughout the rest of the paper. 

A direct measure of the systematic error associated with 
our ED solution is the comparison between right-hand side 
and left-hand side of Eqs. ([8]l. In Ref. 12, using up to Ng — 6, 
the deviation between "lattice" and "impurity" estimates of 
the density, obtained by integrating over the frequency respec- 
tively the right-hand and left-hand sides of the first of Eq. (O 
has been found to strongly increase as the chemical potential 
approaches the bottom of the band and the dilute region is 
approached. Here we found that a careful evaluation of Mat- 
subara sums, together with the different values of /3 and of 
the maximum frequency used in the fit, can make the devi- 
ations much smaller than what found in Ref. 12 already for 
Ns — 6. Increasing Ng to 8 and 9, we have been able to 
reduce the difference between the two estimates down to 1 % 
at n ~ 0.002 — 0.003, virtually eliminating discretization er- 
rors. Therefore the only limitation of our data is the DMFT 
approximation and no further uncertainty is introduced by the 
numerics. 



C. Zero density limit and Universal behavior 

Even if we have almost eliminated the discretization errors, 
a more profound limitation appears in the zero density limit 
when the universal behavior is addressed. The key requisite 
to approach the unitary limit is the divergence of the s-wave 
scattering amplitude Ug- From a diagrammatic point of view 
Ug is obtained as a ladder sum for the irreducible vertex part 
in the s-wave channel at g = in the vacuum. This leads 
to the simple condition for which ag diverges at the smallest 
interaction value for which the two-body problem develops a 
bound state. Within DMFT the only contribution to the ladder 
sum comes from the local vertex functionii. It is straightfor- 
ward to verify that this local approximation is not sufficient to 
recover the divergent Ug, but only leads to a finite value which 
depends on the chosen lattice. However, at static MF level in 
the broken-symmetry superfluid solution a static on-site ver- 
tex function U ^ Uc is sufficient to reproduce the required 
criticality of the two-body problem, despite the non-divergent 
scattering length. As a consequence of the finite "local" scat- 
tering length a'/'^, if we reduce the density maintaining the 
coupling strength at the "universal value" Uc for which the 
full solution of the model would lead to a divergent ag, we 
will eventually enter a regime in which the quantity kpa^"'^ 
decreases, since the vanishing Fermi momentum is not com- 
pensated by the infinite scattering length. This will inhibit to 
reach a real universal regime beyond MF. The necessity of 
non-local contributions in the zero-density limit is even more 
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FIG. 3: (Color Online) Particle density n (solid blue line) and superfluid 
order parameter A (red dashed line) as a function of the chemical potential ^ 
for a semicircular DOS. The inset shows the low density region in logaiithmic 
scale, which shows that n oc fi'^^^ and A oc ^. 

transparent if the symmetry is not broken (normal state). Here 
it is indeed necessary to include a nonlocal vertex function 
(divergent in the ultraviolet in the continuum limit) in order 
to obtain a finite contribution in the fc^ = limit. The above 
analysis implies that, while the DMFT calculations will be ac- 
curate at finite density, when the density becomes extremely 
small the approach will finally acquire a MF character with 
vanishingly small (weak-coupling) corrections. The value of 
the density below which this effect will take place is clearly 
a function of the effective "local" a^, which measures how 
far from the unitary limit we are in practice. This limit den- 
sity can be estimated by simply requiring that kpaJ's'^ ^ 1- 
For the semicircular DOS we obtain a local scattering length 
a^.°= ~ 1.08892. This impHes that DMFT will reduce to static 
MF only for n ^ 1.3 - 10^^, a remarkably small value, smaller 
than the lowest densities attainable in present-time QMC sim- 
ulations. Therefore DMFT is expected to be accurate down to 
extremely low density, making it a reliable tool for the inves- 
tigation of the small-density regime. 

V. DMFT RESULTS 

Before discussing the role of the different lattices in the 
finite-density corrections to the universal limit, in this section 
we analyze the behavior of the most relevant observables in 
our DMFT solutions. For the sake of definiteness, in this sec- 
tion we limit ourselves to the semicircular DOS, which is the 
most easily implemented. 



As we already mentioned, DMFT is naturally expressed in 
the grandcanonical ensemble where the chemical potential is 
the natural variable. In Fig. [3]we plot n and A as a function 
of /i. Starting from half-filling (n = 1), where the value of 
the chemical potential is fixed by the particle-hole symmetry 
of the model Ai(n)|„=i = D + f/c/2 = D- D/2 ^ D/2, the 
density vanishes as the chemical potential vanishes, i.e. ap- 



FIG. 4: (Color Onhne) Ground state energy density (Total energy divided by 
the density) as a function of the density n. We plot the total energy density 
e (blue solid line), the kinetic contribution k (red dashed line), the absolute 
value of the potential energy contribution \v\ (dot-dashed green line). The 
inset shows the low-density region in logarithmic scale. 



proaches the bottom of the band. A first look at the curve 
shows a rather smooth evolution from a high-density region in 
which n has essentially constant slope (i.e., the compressibil- 
ity K = dn/dfi is nearly constant) to a low-density regime 
in which ?i ~ /i"^/^ or k ^ /i^/^ as expected from universal- 
ity. The power-law behavior is emphasized by the logarithmic 
scale used in the inset of Fig. [3] Analogously A evolves from 
being nearly independent on ji for large densities to a linear 
behavior at small /i. 

Another relevant quantity in the unitary limit is the energy. 
To obtain the internal energy E = V + K, we computed the 
potential (V) and kinetic (K) energies. The potential energy 
per site is 

V^ = t/^^KTn,.^) =f/d, (9) 

i 

where N is the total number of lattice sites and d is by defi- 
nition the fraction of doubly occupied sites (local pairs) in the 
ground state, which can be computed directly as a static aver- 
age of the double occupancy operator. The kinetic energy per 
site reads 

keFBZ " 

(10) 

and it can be calculated from the lattice Green's function 
G{e,iujn)- For the semicircular DOS Eq. ( fTOl i is simplified 
and it becomes a function of the local Green's functions only 
(even if a sum over Matsubara frequencies is still required): 

n 

We compared Eq. ( fTol l and (fTTI) obtaining an almost perfect 
agreement, which allows us to confidently use Eq. ( fTOb for 



the other dispersions in the following, for which the simpler 
relation ( fTTT i does not hold. 

In Fig. |4]we plot the energy per particle e — Ejn together 
with the two contributions fc = Kjn and \v\ — |V"|/n. The 
amplitude of the kinetic term is always larger than the poten- 
tial term, correctly implying E ^ K > ^, even if the 
difference rapidly shrinks as the density is decreased. Inter- 
estingly, the two contributions both vanish as v}!"^ for small 
densities, as shown by the logarithmic plot in the inset of Fig. 
m while the total energy, obtained as their sum (or the dif- 
ference between the two curves in Fig. [Hi, scales as r?!'^ , 
the correct behavior for an energy-density in the dilute limit. 
Thesepower-law behaviors have been recently proposed in 
Ref. based on very general arguments, and are confirmed 
by our calculations. To our knowledge, this is the first nu- 
merical confirmation of these results. Moreover, the correct 
behavior of e, which emerges through a subtle cancellation of 
the leading contributions in v and |fc|, is a valuable confirma- 
tion of the accuracy of our estimates for the two contributions 
K and V. 



A. Rescaled Quantities 

In this section we extend the MF analysis of Sec. |III] com- 
puting the same observables using DMFT. As we will detail in 
the following, the dynamical nature of DMFT introduces size- 
able improvements over static MF, but some general trends we 
discussed are found to be robust. 
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FIG. 5: (Color Online) Renormalized chemical potential as a function of 
density. Comparison between static MF and DMFT (marked by open dots). 
Blue dashed line is for the cubic lattice, solid red line for the semicircular 
DOS and dot-dashed green line for the lattice gas. 

The large correction introduced by DMFT with respect to 
MF is evident in Fig. |5l|6]and|7] where we compare, respec- 
tively, the reduced chemical potential y^jEp, superfluid order 
parameter ^jEp and internal energy ^ — EjiijhnEp) ob- 
tained in DMFT and in static MF for the three DOS's. For 
n = 1 the chemical potential is set by particle-hole symmetry, 
and it is therefore the same in DMFT and MF^''. As soon as 
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FIG. 6: (Color Online) RenormaUzed internal energy as a function of den- 
sity. Comparison between static MF and DMFT (marked by open dots). Blue 
dashed Hne is for the cubic lattice, solid red line for the semicircular DOS and 
dot-dashed green line for the lattice gas. 



the density is decreased. Fig. |5]clearly shows that DMFT in- 
troduces a sizeable change with respect to MF, which is due to 
the accurate treatment of quantum fluctuations. Interestingly, 
the improvement brought by DMFT increases by reducing the 
density in a wide range of densities, and it appears much larger 
for the semicircular and lattice gas densities of states than for 
the cubic lattice. The lattice gas displays the weakest depen- 
dence on density. It is also evident, however, that the DMFT 
curves approach the same limit of the static MF as the density 
actually approaches zero, as we expected from the arguments 
given in Sec. II, where we have shown that the including only 
local diagrams the scattering amplitude can not diverge and 
contributions beyond MF are bound to vanish for rt ^ 0. The 
reduction of DMFT to static MF is confirmed by the behavior 
of the self-energies (not shown). While for every finite den- 
sity the self-energies have a non trivial frequency dependence, 
in the very small density range in which DMFT rapidly col- 
lapses on the static MF, this frequency dependence disappears 
and the self-energies become constants, as in static MF. The 
results for the normalized internal energy ^ = Ej {3/5nEp) 
are reported in Fig. |6] and show similar trends with respect 
to the chemical potential, confirming the anomalous behavior 
of the cubic lattice and the weaker density dependence of the 
lattice gas. 

In Fig. |2]we propose a similar comparison also for the su- 
perfluid gap A, again divided by the natural energy scale Ep. 
As we have akeady discussed for the MF results, here the 
three different lattices have a similar behavior All the curves 
are monotonically decreasing functions of the density, in con- 
trast with the minima presented by the chemical potential and 
the energy in some of the lattices. However, DMFT deter- 
mines strong corrections also for this quantity. It is intriguing 
that the different lattices receive significantly different renor- 
malizations. The cubic lattice, that has the largest size effects 
in static MF is the one that is less affected by the dynamical ef- 
fects introduced by DMFT, while the semicircular and lattice 
gas DOS have similar shifts. In all cases the superfluid gap 
is reduced, as expected by the introduction of effects beyond 
MF. Nonetheless, also for A the zero density limit becomes 
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FIG. 7: (Color Online) Renormalized superfluid order parameter as a func- 
tion of density. Comparison between static MF and DMFT (marked by open 
dots). Blue dashed line is for the cubic lattice, solid red line for the semicir- 
cular DOS and dot-dashed green line for the lattice gas. 



identical to the static MF. 

Nonetheless, it is indeed evident that, for all the observ- 
ables, the improvement brought by DMFT on MF increases 
as the density is reduced, until we reach a very low density for 
which the non-divergence of the scattering amplitude forces 
the result to rapidly approach the static MF. In this light, it is 
safe (and even prudentic) to assume that down to the densities 
at which the distance between MF and DMFT data is still in- 
creasing in Figs l5l6l'7l DMFT data are not plagued by the lack 
of a divergent as- In the semicircular DOS and the lattice gas 
and for ji/Ep and ^ this density scale is around n ~ 0.02. 

A clear outcome of our DMFT analysis is that the non- 
universal finite-density effects can be very large and they 
strongly depend on the choice of the lattice. For all lat- 
tices, the evolution from the large-density regime to the dilute 
regime is very smooth and regular The difference between 
the different lattices is still very large at n ~ 0.05 0.1 or 
smaller, which are therefore not representative of a real low- 
density regime where the physics becomes universal. This 
implies that any attempt to obtain informations about the uni- 
versal regime should be based on a careful finite-density limit, 
and that densities much smaller than n ~ 0.05-^0.1 should be 
used in this extrapolation. Unfortunately modern QMC calcu- 
lations are just at the limit of this region (To our knowledge 
the smallest density used for the attractive Hubbard model is 
n — 0.05 in Ref. |9|), suggesting that further work is required 
to extract in a completely reliable way the properties of the 
universal Fermi gas. 

A second important result is that the cubic lattice is the 
worst choice among the models we considered. In this model 
the large-density regime has a peculiar behavior which results 
from the singularities of the DOS, and it is very hard to wash 
out these effects by reducing the density. As a result, this 
model has the slowest convergence to the zero-density limit. 
Furthermore, the corrections introduced by DMFT are smaller 
than for the other lattices. 

The lattice gas and the semicircular DOS give rise to com- 
parably accurate results, even if the lattice gas presents a 
weaker dependence on the density, which makes it the best 



FIG. 8: (Color Online) Dependence on the density of the chemical potential 
divided by the lattice free Fermi-energy for the three lattices under study 



candidate for future studies. The correction introduced by 
DMFT on static MF is similar for the two models. 

A direct comparison between our DMFT results and QMC 
simulations is not straightforward in light of the limits of 
DMFT in the zero-density limit. Nonetheless, if we simply 
consider the values of the observables around the density be- 
low which DMFT rapidly approaches the static MF, we find 
that the results of the lattice gas are in very good agreement 
with other estimates. For example for the lattice gas we have 
^ ~ 0.43 and A ~ 0.55 (the estimate for A is rather arbitrary, 
because this quantity does not show any minimum as a func- 
tion of density). Comparing the latter two quantities with Fig. 
14 of Ref. where several results from different numerical 
methods are collected, our data are clearly in the same range. 

We conclude our analysis be presenting, as an example, the 
DMFT data for the chemical potential normalized by the lat- 
tice Fermi energy corresponding to each DOS. This shows that 
also the cubic lattice displays now a minimum as a function 
of density, which was not present in Fig. |5] Indeed the com- 
parison with the free-space Fermi energy is only meaningful 
at small density, especially for the cubic lattice, whose DOS 
is very different from the free one (see Fig. [T]). This plot 
shows that the existence of a minimum is a general results, 
which can be used to define a safe region in which DMFT is 
an improvement on static MF and an accurate description of 
a finite-density lattice model. This makes the argument about 
"fake universality" that we discussed in Sec. |lll] relevant for 
any finite-density extrapolation to the zero-density limit. 



VI. SUMMARY AND CONCLUSIONS 

We have presented a thorough investigation of the approach 
to the zero-density limit of a lattice model with local attrac- 
tive interaction. Three different lattice models have been con- 
sidered, all sharing the low-energy density of states of free 
fermions in three dimensions. Tuning the interaction strength 
at the unitary limit, for which the scattering length diverges. 
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the three models approach the same zero-density limit, which 
is expected to reproduce the properties of a three-dimensional 
unitary Fermi gas. 

Our investigation is based on Dynamical Mean-Field The- 
ory (DMFT), and for comparison on static mean-field at finite 
density. We show that DMFT introduces large corrections to 
MF down to very small densities, even if it is bound to reduce 
to static MF at zero density because the local approximation 
inherent to DMFT does not allow for a divergence of the scat- 
tering length beyond MF. Following the evolution from large 
densities of the order of one fermion per site to zero density, 
the improvement introduced by DMFT increases as the den- 
sity is reduced for a large window of densities. Only at very 
small densities of the order of 10^^ DMFT starts to approach 
static MF. The finite-density corrections are non-universal as 
expected and they turn out to strongly depend on the choice 
of the lattice. For all the lattice we considered, the evolu- 
tion as a function of density is really smooth, and the large- 
density properties influence the physics down to very low den- 
sities. The three-dimensional cubic lattice, in particular, has 
the slowest convergence because of its pecuhar DOS. The lat- 
tice gas and a system with a semicircular DOS display a much 



better convergence. 

While the chemical potential and the energy strongly de- 
pend on the chosen lattice, the superfluid order parameter has 
a more lattice-independent behavior, even if the finite-density 
effects are sizeable. This can be associated to the fact that the 
unitary regime is closer to a strong-coupling Bose-Einstein 
regime, in which the interaction strength controls the order 
parameter and the details of the bandstructure are less impor- 
tant. 

Our result highhght the importance of non-universal finite- 
density corrections to the unitary limit and impose serious 
constraint on any attempt to extrapolate the physics of the uni- 
versal Fermi gas from finite-density calculations. In particular 
densities around n ~ 0.05 0.01 are not representative of the 
dilute regime, but they are still controlled by the large-density 
physics (and they strongly depend on the chosen lattice). 

We finally remind that the limitations of DMFT are only 
dangerous when we want to describe the zero-density Umit of 
a f ermionic gas at unitarity. When we directly deal with lattice 
models, our results shown that the approach proves extremely 
accurate and introduces large corrections to static mean-field 
down to extremely low densities. 
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